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these interactions are determined from the empirical hyperon single-particle potentials in symmetric 
nuclear matter at saturation density. The extended MDI interaction is then used to study in the 
mean-field approximation the equation of state of hypernuclear matter and also the properties of 
hybrid stars by including the phase transition from the hypernuclear matter to the quark matter 
at high densities. In particular, the effects of attractive and repulsive EN interactions and different 
values of symmetry energies on the hybrid star properties are investigated. 
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I. INTRODUCTION 

The study of the in-medium baryon-baryon efFective 
interactions is one of the fundamental problems in nu- 
clear physics. While significant progress has been made 
in understanding the nucleon-nucleon interaction in nu- 
clear medium, our knowledge on the hyperon-nucleon and 
hyperon-hyperon interactions in nuclear medium is very 
limited. The latter is important for understanding the 
properties of hypernuclei, the equation of state (EOS) of 
dense baryonic matter, the properties of compact stars in 
which hyperons are expected to appear abundantly, and 
strangeness production in heavy-ion collisions. There- 
fore, it is of great interest to develop an effective model 
for the hyperon-nucleon and hyperon-hyperon interac- 
tions in nuclear medium. 

The in-medium single-particle potential of a nueleon 
generally depends not only on the nuclear density and 
its momentum but also on its isospin. The isospin- and 
momentum-dependent MDI interaction is such an effec- 
tive nuclear interaction based on the finite-range Gogny 
interaction This interaction has recently been ex- 
tensively used in studying the isospin and momentum- 
dependent effects in nuclear reactions and theproperties 
of neutron stars (for a recent review see Ref. [3 ) ■ In par- 
ticular, using this interaction in the Isospin-dependent 
Boltzmann-Uheling-Ulenback (IBUU) transport model 
with the isospin-dependent in-medium nucleon-nucleon 
scattering cross sections, a relatively stringent constraint 
on the density dependence of the nuclear symmetry en- 
ergy at subsaturation densities was obtained from the 
isospin diffusion data in intermediate-energy heavy-ion 
collisions 0"|5| . The resulting symmetry energy has been 
used to constrain the neutron skin thickness of heavy nu- 
clei @ and the properties of neutron stars 04111 includ- 



ing the transition density, which separates the liquid core 
from the inner crust of a neutron star, and other proper- 
ties of neutron stars by assuming that they only consist 
of nucleons, electrons, and muons [l^. [l3|. 

In this paper, we extend the MDI interaction to include 
nucleon-hyperon and hyperon-hyperon interactions. This 
is achieved by assuming that the nucleon-hyperon and 
hyperon-hyperon interactions have the same density and 
momentum dependence as the nucleon-nucleon interac- 
tion with the interaction parameters fitted to known ex- 
perimental data at normal nuclear matter density. With 
this extended MDI interaction, we then study, as an ex- 
ample, the properties of hybrid stars that not only con- 
sist of appreciable fraction of hyperons but also a pos- 
sible quark matter in their dense core [l4| - [26| . For the 
hadron-quark phase transition, we use the Gibbs con- 
struction [23, with the quark phase described by a 
simple MIT bag model [igI. [29|. In particular, we study 
the effects of attractive and repulsive EN interactions 
and different symmetry energies on the properties of hy- 
brid stars. We note that there have been extensive stud- 
ies on hybrid stars based on various approaches, such 
as the relativistic mean field (RMF) model pjlMMj . the 
Brueckner-Hartree-Fock theroy |35l - l39l | , and phenomeno- 
logical models (40ll4l|. 

This paper is organized as follows. In Sec. |lTl we de- 
scribe the method used in extending the MDI interaction 
to include the hyperon-nucleon and hyperon-hyperon in- 
teractions. The extended MDI interaction is then used 
in Sec. IIII Al to study the thermodynamical properties 
of hypernuclear matter and its equilibrium conditions. 
The MIT bag model is used in Sec. IIII 5] to describe the 
properties of the quark matter and in See. IIII CI to study 
the hadron-quark phase transition in dense matter. In 
Sec. IIVI we show and discuss the results on the parti- 
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cle fractions in dense matter and its equation of state as 
well as the mass-radius relation of hybrid stars in differ- 
ent scenarios. A summary is given in Sec. |Vl 



II. THE MDI INTERACTION WITH 
HYPERONS 

The MDI interaction is an effective nuclear interac- 
tion with its density and momentum dependence deduced 
from the phenomenological finite-range Gogny interac- 
tion In the mean- field approximation, the poten- 
tial energy density in a nuclear matter of density p and 
isospin asymmetry S = (p„ — Pp)/p with pn and pp being, 
respectively, the neutron and proton densities, is given by 



Vip,d) 



+ (Pn+Pp) + TT^" 

Po 2po a + 1 pI 



(1) 



where t{t') is the nucleon isospin, /r(?^,p) is the nucleon 
phase distribution function, and po = 0.16 fm^'^ is the 
saturation density of normal nuclear matter. Values of 
the parameters Ai{x), B, a. A, Ci = Ct,t and 

Cu = Ct.-t can be found in Refs. [ll, The parameter 
X, related to the coefficient of the spin-exchange operator 
in the Gogny/Skyrme interactions, is used to adjust the 
density dependence of the symmetry energy away from 
the saturation density without changing the properties of 
symmetric nuclear matter. At subsaturation densities, its 
value has been constrained between and —1 from the 
analysis of the isospin-diffusion and neutron skin thick- 
ness data i-i,!!!. 



A. The extended MDI interaction 

To extend the MDI interaction for the nucleon-nucleon 
[NN) interaction to include the nucleon-hyperon {NY) 
and hyperon-hyperon (YY) interactions, we assume that 
the latter have the same density and momentum depen- 
dence as that between two nucleons. The potential en- 
ergy density of a hypernuclear matter due to interactions 
between any two baryons then has the following general 
form: 
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where b (&') denotes the baryon octet included in the 
present study, i.e., N, A, S, and S. We use the conven- 
tions that r^v = — 1 for neutron and 1 for proton [43l |. 
TA = for A, TE = -1 for S", for and 1 for S+, 
and Ts = — 1 for 5~ and 1 for S". In the above, the to- 
tal baryon density is now given by p = J2b X^r^ P^b > ^^"^ 
frbi^iP) is the phase-space distribution function of par- 
ticle species t^. The interaction parameters are denoted 
by Ahb' ; ^ftf,'' ^bb', and Cr^^r^- If there are only nucleons, 
we can rewrite Ann = {Ai+Au)/2, A'^^^ = (A;-^„)/2, 
Bnn = B, and Cr^^r'„ = Q for tn = and = C„ for 
tn ^ t'j^, which then reduce to the original parameters 
in the MDI interaction for nucleons The param- 

eter X is again used to model the isospin effect on the 
interaction between two baryons, and its value is taken 
to be or —1 in the present study for all baryon pairs. 

The single-particle potential for a baryon of species Tf, 
in a hypernuclear matter can then be obtained from the 
total potential energy density of the hypernuclear matter, 
given by Vhp = (1/2) Y.b b' ^bb' , as 
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For the interaction parameters Ai,b' , A^i,/ , -Bfet' and 



that involve hyperons, they can in principle 



be determined from the nucleon-hyperon and hyperon- 
hyperon interactions in free space. Because of the 
lack of hyperon-nucleon scattering experiments, infor- 
mation on the hyperon-nucleon interactions has been 
mainly obtained from the hyperon single-particle po- 
tentials extracted empirically from studying A [i^ as 
well as E [ill, H^l and S [43l production in nuclear 
reactions. Although this has led to extensive studies 
of the hyperon-nucleon interaction in the past (48l - [50| . 
the isospin and momentum dependence of the hyperon- 
nucleon in-medium interactions is still not very well de- 
termined, and the situation is even worse for hyperon- 
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hyperon interactions. We therefore assume in the present 
study that A^b'^ ^h^, Bbb' and CTb,T'^, are all proportional 
to corresponding ones in the nucleon-nucleon interaction. 
Specifically, for Abb' , A'^^b' ^"^^ Bbb' , we assume 

Abb' = fbb'ANN, 
^'bb' ^ fbb'A'ffj^, 

Bbb' = fbb'BjyN, (4) 

and for C,-. -r' we assume 

f fbb'^^ (Tborr^, =0), 

Cr,.r' = <^ fbb'Cl (Tb=r^, ^0), 

[fbb'Cu (Tb^r^, 7^0), 

with A and treated differently. 

We determine the values of fbb' by fitting the empirical 
potential uj^'^ ^ of baryon b at rest in a medium consisting 
of baryon species b' . For hyperons in symmetric nuclear 
matter at saturation density, their potentials are 

t^r^(PAr = Po)--30MeV (5) 

for the A potential from the analysis of (7r+,i^+) and 
(A'~,7r~) reactions on nuclei [Uls^l and 

4'^^W = Po) = -18MeV (6) 

for the S potential from the analysis of (S, \H) [H^ and 
{K^,K~^) [s^, HE] reactions. This leads to the values 
/tva = 1 and Jns — 0.73. For the E hyperon, its po- 
tential was taken to be attractive in earlier studies [45|, 
but more recent analysis indicate that it is repulsive in 
the nuclear medium f56l - l60j . In the present work, we 
consider both the attractive and repulsive cases 

U^''\pn = Po) = ±30MeV. (7) 

By setting /jvs = 1 wc get an attractive SN interac- 
tion, and this is called MDI-Hyp-A in the following. To 
obtain a repulsive SN interaction, called MDI-Hyp-R 
in the following, we adjust the values of positive and 
negative terms in the single-particle potential by setting 
Bjvs = 2.33i3jv7v and A'j^^: = 2.33A^jy without chang- 
ing other parameters as in the case of MDI-Hyp-A in- 
teraction. A similar method of changing an attractive 
SN interaction to a repulsive one was used in the RMF 
calculation (sij by changing the coupling constants of lo 
and p mesons. For the hyperon-hyperon interaction, the 
parameters are fitted according to [6l| 

=Po)^-40MeV, (8) 

which gives the strength of the hyperon-hyperon interac- 
tions as /aa = 0.63, /as = 1-25, /a3 = 1.25, /ss = 0.49, 
/sH = 1, and /hh = 0.53. Detailed values of the param- 
eters are listed in Table. HI These parametcrizations can 
be viewed as a baseline for studying the properties of the 
hypernuclear matter, and more sophisticated treatments 
are left for future work after the in-medium properties 
of hyperons are better understood. It will be shown in 
the following that many interesting results can already 
be obtained even with these simple parametcrizations. 



B. Single-particle potentials 

An important quantity related to the interaction of a 
particle in nuclear medium is its single-particle potential 
as given by Eq. which is also needed later in our study 
of the properties of neutron stars. The single-particle po- 
tential of a particle depends not only on the density of the 
medium but also on the momentum of the particle. In 
this section, we show and discuss the single-particle po- 
tentials of both nucleons and hyperons in nuclear matter 
obtained from the extended MDI interaction. 



MDI-Hyp interaction ' 
for symmetric / 
nuclear matter ' 




FIG. 1: (Color online) Single-particle potentials for particles 
at rest in symmetric nuclear matter as functions of density. 

We first show in Fig. [1] the single-particle potential of a 
particle at rest in symmetric nuclear matter as a function 
of density. Although the nucleon potential is more attrac- 
tive at normal density than those of hyperons, it becomes 
more repulsive than the hyperon potentials above about 
5 times normal nuclear density, including the S potential 
that is attractive at normal density. For the S poten- 
tial that is repulsive at normal density, it becomes more 
repulsive as the density increases and becomes slightly 
attractive only at very low densities. Comparing our 
results with those from other models given in Ref. [g^I 
(and references therein), the single-particle potentials of 
A and E are close to those from the chiral effective field 
theory [g^l, but more repulsive than those based on the 
G- matrix calculations using the soft core Nijmegen model 
or the Jiilich meson-exchange model for the free hyperon- 
nuclcon interactions, particularly at high densities. 

The single-particle potential of a particle from the ex- 
tended MDI interaction also depends on its momentum. 
Figure O shows the momentum dependence of the single- 
particle potentials for both nucleons and hyperons in 
symmetric nuclear matter at saturation density. Again, 
results using both attractive and repulsive SN interac- 
tions (i.e., MDI-Hyp-A and MDI-Hyp-R) are shown for 
comparison. Also indicated in the figure is the Fermi 
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TABLE I: Parameters for the MDI-Hyp-A and MDI-Hyp-R interactions with a; = and a; = — 1. All except a are in units of 
MeV. A'j^Y,{A.) and Bmi:{A) are for the MDI-Hyp-A interaction, and A'j^^{R) and Bns{R) are for the MDI-Hyp-R interaction. 
Other parameters are the same for both interactions. 
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FIG. 2: (Color online) Single-particle potentials in symmetric 
nuclear matter at saturation density po as functions of particle 
momentum. The Schrodinger equivalent potential obtained 
by Hama et al. [6^ . [6^ from the nucleon-nucleus scattering 
data is shown by stars for comparison. E(A) and S(R) are for 
the MDI-Hyp-A and MDI-Hyp-R interactions, respectively. 



momentum of nucleons. For nucleons, the single-particle 
potential from the MDI interaction is consistent vifith the 
Schrodinger equivalent potential obtained by Hama et al. 
from the nucleon-nucleus scattering data jll, HH up to 
the nucleon momentum of 500 McV. For hyperons, the 
momentum dependence of their single-particle potentials 
from the extended MDI interaction is similar to that ob- 
tained from the G-matrix calculations based on the free 
Nijmegen NY interaction [STi]. Both show an increase 
with increasing momentum and are similar at low mo- 
menta as both are constrained by available experimen- 
tal data. They are, however, slightly different at high 
momenta. The momentum dependence of NY and YY 
interactions thus remains an open question, especially 



at high momenta. Also, the density dependence of the 
high-momentum behavior of the mean-field potential is 
poorly known. Since the maximum Fermi momenta of 
nucleons and hyperons reached in hybrid stars are not 
very high, the incorrect nucleon potential and the un- 
certainty of hyperon potentials at high momenta given 
by the extended MDI interaction are thus not expected 
to affect significantly the properties of cold hypernuclear 
matter considered in the present study. 
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FIG. 3: (Color online) Symmetry potentials of nucleon, E 
and H in asymmetric nuclear matter at saturation density 
p = po. T.{A) and E(7?) are for the MDI-Hyp-A and MDI- 
Hyp-R interactions, respectively. 

It is known that both proton and neutron single- 
particle potentials in asymmetric nuclear matter of un- 
equal proton and neutron densities are approximately lin- 
ear in the isospin asymmetry 5 of the matter. Whether 
this is also the case for hyperons is not clear in the litera- 
ture. For the extended MDI interaction introduced in the 
present study, such a linear dependence in isospin asym- 
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metry, however, also holds for E and S hyperons. The 
single-particle potential of a particle in asymmetric nu- 
clear matter can thus be written in general as Ur,,{p, S) w 
Ur^{p,S = 0) — TbU^y^{p)S in terms of the symmetry 
potential U^ynAp)^ defined by U^„^{p) = {Un{p,S) - 



Up{p,5))/2S, U^yM 



{U^-ip,S) -Uj:+{p,6))/26 and 



C^sym(p) = {Ue-{p,S) ~ C/ho(/3,(5))/2(5 for the nucleon, S 
and S hyperons, respectively. From the single-particle 
potentials of nuclcons, S and S hyperons in asymmet- 
ric nuclear matter at normal nuclear matter density and 
of isospin asymmetry S = 0.2, we have calculated their 
symmetry potentials using above definitions. In Fig. [31 
the momentum dependence of these symmetry potentials 
are compared. All symmetry potentials are seen to de- 
crease with increasing momentum. At zero momentum, 
the symmetry potentials of nucleon and S arc about 34 
MeV and 12 MeV, respectively, and for the E hypcron 
they are 17 MeV for an attractive EN interaction and 9 
MeV for a repulsive one. 
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FIG. 4: (Color online) Symmetry potentials of nucleon, E and 
H in asymmetric nuclear matter at density p = 3pQ. and 
E(7?) are for the MDI-Hyp-A and MDI-Hyp-R interactions, 
respectively. Note that different scales are used for a:: = (left 
panel) and x = —1 (right panel). 

Although the symmetry potentials at normal nuclear 
density are independent of the value of x, which is used 
in the MDI interaction to model the stiffness of nuclear 
symmetry energy at densities different from the normal 
density, this is not the case at other densities. This is 
demonstrated in Fig. |4] for the symmetry potentials of 
the nucleon, the E hyperon and the S hyperon in asym- 
metric nuclear matter at density p = 3po for the two sym- 
metry energy parameters a; = and x = —1. For x = 0, 
the symmetry potentials at zero momentum are 53 MeV 
and 19 MeV for the nucleon and the S hyperon, respec- 
tively, and for the E hyperon they are 26 MeV with the 
MDFHyp-A interaction and 2 MeV with the MDFHyp-R 



interaction. These values are changed to 174 MeV and 63 
MeV for the nucleon and the 2 hyperon, respectively, and 
for the E hyperon to 87 MeV for the MDI-Hyp-A inter- 
action and 143 MeV for the MDI-Hyp-R interaction for 
the case of a; = — 1 . As noted in Ref . [66[ , the charged E 
baryon ratio in heavy-ion collisions can be used as a probe 
to constrain the nuclear symmetry energy (potential) at 
densities larger than 3po- It will be very interesting to see 
how the symmetry potentials of S and S hyperons in nu- 
clear matter affect the charged E hyperon ratio and the 
charged S hyperon ratio in heavy-ion collisions induced 
by neutron-rich nuclei. This is important for determin- 
ing the high density behavior of the nuclear symmetry 
energy using these ratios in heavy-ion collisions. The ex- 
tended MDI interaction with hyperons is therefore useful 
for studying the nuclear symmetry energy (potential) at 
p > 'ipo in the transport model calculation for heavy- ion 
collisions. 



III. EQUILIBRIUM CONDITIONS AND 
THERMODYNAMICAL RELATIONS IN DENSE 
MATTER 

The extended MDI interaction is also useful for study- 
ing the properties of a hybrid star that is expected to 
have a quark core at high densities, a mixed phase of 
quarks and hadrons at moderate densities, and a hadron 
phase at low densities. We review in this section the /3- 
equilibrium, charge neutrality, and baryon number con- 
servation conditions of such a matter. 



A. The hadron phase 

In the hadron phase of a hybrid star, the matter con- 
sists of nucleons, hyperons and leptons. For leptons, we 
include both electrons and muons with their masses taken 
to be and 106 MeV, respectively. At equilibrium, these 
particles satisfy following baryon number conservation, 
charge neutrality and /3-equilibrium conditions: 



^Pibi = P, 

i 

^Ptli = 0, 

i 



(9) 
(10) 
(11) 



In the above, qt and bi are, respectively, the charge and 
baryon numbers of particle species i, where i can be nu- 
cleons, hyperons or leptons, and their density and chemi- 
cal potential arc denoted, respectively, by pi and pi. The 
total baryon density of the hadron phase is denoted by p, 
and p^ and p^ are, respectively, the baryon and charge 
chemical potentials of the hadron phase. 

Taking into account their interactions in the mean-field 
approximation, the chemical potential of baryon species 
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Tf, is given by 



2m, 



(12) 



where Pp^^ is their Fermi momentum, and Ur,, and m-r^ 
are, respectively, the single-particle potential and mass 
of the baryon. For leptons {I = fj., e) their chemical 
potential is given by fii = {mf + p^iY^^ with = 
(Stt^P;)^/'^ being their Fermi momentum. The relative 
abundances of various hadrons and leptons for a given 
total baryon density are then obtained by solving above 
equations. 

In terms of the densities of various particles, the total 
energy density of the hadron phase can be written as 



Vh + Vl, 



(13) 



where Vh and Vl are the contributions from baryons and 
leptons, respectively. The former can be written as 



Vh — Vhp + Vhk + Vhm , 



(14) 



where Vhp = (1/2) J2b w ^bb' is the potential energy den- 
sity of baryons with Vw calculated from Eq. ([2]), and 
Vhk and Vhm are, respectively, the kinetic energy and 
mass contributions given by 



Vhk = iOttC 

b Tb 
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(15) 
(16) 



The contribution Vl from leptons is calculated by treat- 
ing them as a free Fermi gas, i.e.. 
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The pressure of the hadron phase is obtained from the 
thermodynamical relation 

= + = (^ ^ - Vh) + (Y, pipi - Vl). 

b ; 

(18) 

where b and / run over all species of baryons and lep- 
tons, respectively. We note that the thermodynamical 
consistency condition 
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dp 



(19) 



B. The quark phase 

As the nuclear matter density increases, such as in the 
core of neutron stars, not only hyperons appear but also a 
quark matter could exist [l5| . To take into consideration 
possible transition between the hadronic matter and the 
quark matter, we follow many previous studies to use in 
the present study the MIT bag model [IB, [H] to describe 
the cold quark matter that might exist in the dense core 
of neutron stars. 

For the quark phase consisting of quarks and leptons, 
the baryon number conservation and charge neutrality 
conditions are given by expressions similar to Eqs. 
and (jlOp with i denoting now quarks and leptons. For 
the /3-equilibrium condition in the quark phase, it is given 
by pLi = p^hi — n^Qi with p'^ and p^ being the baryon 
and charge chemical potentials of the quark phase, re- 
spectively. 

The total energy density and pressure of the quark 
phase can be calculated from 

e'^ = VQ + VL, (20) 

P« = Pq + Pl, (21) 

where Vq and Pq are the energy density and pressure 
of quarks, which can be calculated from the MIT bag 
model as described in the following, and Vl and Pl are 
the energy density and pressure of leptons given by same 
expressions as those in the hadron phase. 

At zero temperature, the density pf, chemical potential 
/iy, and energy density Vf of quarks of flavor f — u,d, s 
in the quark matter are given, respectively, by 
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(22) 



where ppj^ is the Fermi momentum of quarks of flavor /. 
For the quark masses, they are taken to be m„ = to^ = 
and TTis = 150 MeV. In the bag model, the energy density 
is modified by a bag constant B, resulting in an energy 
density given by 



Vq 



f 



B. 



(23) 



This leads to the following pressure for the quark matter 

id{VQ/p) 



Pc 



f 



Vq=p' 



dp 



(24) 



where p is the total baryon density of the quark phase 

lj2pf- (25) 



P = 



is satisfied. 
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C. The hadron-quark phase transition 

The hadron-quark phase transition leads to a mixed 
phase of hadronic and quark matters, which is usually 
described by the Gibbs conditions [23, [28j 

rpH _ rjiQ pH _ pQ 

fJ'b = f^b = M?, Mc=Mf = M?- (26) 

The Gibbs conditions for the chemical potentials can also 
be expressed as 

_ 1 2 

= A^s = ^A*n + ^A^e- (27) 

Since only the case of zero temperature is considered in 
this paper, the first condition in Eq. ([26|) is always satis- 
fied. To solve above equations, we follow the method of 
Ref. [2^. In this method, one first calculates the pressure 
of the hadron phase at a series of baryon densities. From 
the known chemical potentials fi„ and fi^ of the hadronic 
matter at these densities, one then calculates the pres- 
sure of the quark matter that is in chemical equilibrium 
with corresponding hadronic matter, i.e., with the quark 
chemical potentials fiu, fJ-d and fig determined by the two 
relations in Eq. (|27p . The hadronic baryon density at 
which the two pressures are the same is then the low- 
density boundary of the hadron-quark phase transition. 
As an illustration, we show in Fig. [5] the results for the 
case of the symmetry energy parameter a; = and the 
bag constant B^/^ = 180 MeV. The solid line is the 
density dependence of the pressure of the hadron phase 
calculated using the model described in Sec. IIII A| while 
the dash-dotted line is the pressure of the quark matter 
that is in chemical equilibrium with the hadronic mat- 
ter as a function of the baryon density of the hadronic 
matter. The low-density boundary of the hadron-quark 
phase transition is indicated by the intersection of the 
solid line and the dash-dotted line. 

Above the low-density boundary determined above, 
the dense matter enters the mixed phase, in which the 
hadron phase and the quark phase satisfy following chem- 
ical equilibrium, mechanical equilibrium, baryon number 
conservation, and charge neutrality conditions: 

^ibh ~ iJLcqi = Mi, 

pH ^ pQ^ 

b rt f 

(1 - ^) E E + 1 E pfif + E = 0' 

(28) 

where i runs over baryons. leptons and quarks, and Y 
is the baryon number fraction of the quark phase. We 
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FIG. 5: (Color online) Pressures as functions of baryon den- 
sity p for the hadron phase (solid line), mixed phase (dotted 
line) and quark phase (dashed line). The pressure of the quark 
matter that is in chemical equilibrium with the hadronic mat- 
ter as a function of the baryon density of the hadronic matter 
is denoted by the dash-dotted line, whose intersection with 
the solid line gives the low-density boundary of the hadron- 
quark phase transition. 

note that in the mixed phase the total charge is zero 
and the leptons play an important role in maintaining 
the charge neutrality and /3-equilibrium conditions. The 
total energy density and pressure of the mixed phase are 
calculated according to 

- il-Y)VH + YVQ + VL, 
= {1-Y)Ph + YPq + Pl. 

(29) 

It is obvious that we have a pure hadron phase for F = 
and a pure quark phase for F = 1. For the case of the 
symmetry energy parameter a; = and the bag constant 
B'^/'^ = 180 MeV considered above, the pressure of the 
mixed phase as a function of the baryon density obtained 
from solving Eq. ((28|) is shown in Fig.[5]by the dotted line. 
The mixed phase starts at the low-density boundary of 
the hadron-quark phase transition and ends at the high- 
density boundary of the hadron-quark phase transition 
when the matter is a pure quark matter, whose pressure 
is calculated using the model described in Sec. IIII Bl and 
is shown by the dashed line in Fig. [S] as a function of the 
baryon density. 

There is a recent study based on the Relativistic Mean- 
Field model for the hadron phase and the MIT bag model 
for the quark phase to compare the behavior of the mixed 
phase in the Gibbs construction with that in the Maxwell 
construction, which does not require same charge chem- 
ical potential for the hadronic matter and quark matter 
in the mixed phase (67| . The pressure of the mixed phase 
in the Maxwell construction is found to be constant with 
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respect to its baryon density, which is in contrast with 
the increasing pressure of the mixed phase with increas- 
ing baryon density in the Gibbs construction as seen in 
Fig. [5] It is worthwhile to point out that although the 
charge neutrality condition in Eq. ((28| is satisfied glob- 
ally, it is violated in each phase. A more realistic equa- 
tion of state of the mixed phase can be obtained from 
the Wigncr-Scitz cell calculation by taking into account 
the Coulomb and surface effects. The equation of state 
of the mixed phase obtained from this approach lies be- 
tween those from the Gibbs and Maxwell constructions. 
The latter can thus be viewed as two extreme cases cor- 
responding to certain values of the surface tension pa- 
rameter in the Wigner-Seitz cell calculation [il, |6§| . 



IV. RESULTS AND DISCUSSIONS 

In this section, we use the extended MDI (MDI-Hyp) 
interaction to study the equation of state and the rela- 
tive particle fractions in the charge neutral and /3-stable 
hypernuclear matter. Including the hadron-quark phase 
transition in the hypernuclear matter, we further study 
the mass-radius relation of hybrid stars. Results from 
different values of the symmetry energy parameter a;, the 
bag constant, and attractive and repulsive SN interac- 
tions are compared. 



A. The hypernuclear matter 
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FIG. 6: (Color online) Particle fractions in a hypernuclear 
matter for the MDI-Hyp interaction with a; = ((a) and (c)) 
and X — —1 ((b) and (d)) as well as with attractive ((a) and 
(b)) and repulsive ((c) and (d)) EN interactions. 

Using the extended MDI interaction with hypcrons 
(MDI-Hyp), we have studied the particle fractions in a 
hypernuclear matter as functions of the baryon density. 




0246802468 10 

P/Po 

FIG. 7: (Color online) Same as Fig.[6]for the chemical poten- 
tials in a hypernuclear matter. 



The results are shown in Fig.[n]for MDI-Hyp-A ((a) with 
X = and (b) with .t = -1) and MDI-Hyp-R ((c) with 
X = and (d) with x — —1) interactions. The density 
dependence of corresponding particle chemical potentials 
is shown in Fig. [71 It is seen that the effect of hypcrons 
is more pronounced at higher densities, and the appear- 
ance of hypcrons prevents the neutron chemical potential 
from increasing too fast. Similar to the results from other 
works, the A hyperon appears at baryon density of about 
3po- For the S hyperon, the density at which it appears 
depends on the sign of the SN interaction. For the at- 
tractive MDI-Hyp-A interaction, it appears at about 2po, 
while for the repulsive MDI-Hyp-R interaction, it does 
not appear until very high densities due to the rapid 
increase of the chemical potentials of S hypcrons with 
increasing total baryon density as shown in Fig. [3 mak- 
ing it difficult to satisfies the /3-equilibrium condition at 
small baryon densities. This result is similar to that re- 
ported in the literature [Sll, H^. Also, with the value 
X — corresponding to a stiffer symmetry energy, the 
fractions of Icptons are lower at subsaturation densities 
but higher at higher densities, which results in a larger 
charge chemical potential at higher densities and the ap- 
pearance of negatively (positively) charged hypcrons at a 
lower (higher) density. A stiffer symmetry energy also in- 
creases the neutron chemical potential and thus leads to 
a larger total hyperon fractions at higher densities. Since 
the total charge of the hypernuclear matter is conserved, 
an increase in the fraction of negatively charged hyperon 
results in an increasing fraction of protons and decreas- 
ing fraction of leptons. We note that for the MDI-Hyp-R 
interaction and at the density lOpo, about 50% of the 
particle fraction is made up by hyperons, and the frac- 
tion of A particles is larger than that of protons at high 
densities. 

In Fig. [HI we show the equation of state of nuclear mat- 
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B. The hadron-quark phase transition 
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FIG. 8: (Color online) Equation of state of pure nuclear mat- 
ter and hypemuclear matter from the MDI-Hyp interaction 
with a; = and s = — 1 as well as with attractive (upper 
panel) and repulsive (lower panel) EN interactions. 



tcr from the MDI interaction with a; = and x — ~\ and 
of hypemuclear matter from the MDI-Hyp-A and MDI- 
Hyp-R interactions. It is seen that the EOS is softened 
when hyperons are present as compared to the EOSs of 
pure nuclear matter. We note that the EOS at mod- 
erate densities plays an important role in determining 
the maximum mass of neutron stars. Without hyperons, 
the EOS at moderate densities is stiffer for smaller value 
of while it becomes stiffer for larger value of x when 
hyperons are included, especially for a repulsive SN in- 
teraction. As shown in Fig. [6l a stiffer symmetry energy 
leads to a larger number of hyperons, which results in a 
softer EOS as a result of the lower pressure due to the 
presence of more degrees of freedom, and the effect of hy- 
perons on the EOS is smaller for a soft symmetry energy, 
as fewer hyperons are then present in the hypemuclear 
matter. One the other hand, the symmetry energy con- 
tribution to the pressure of the hypemuclear matter is 
larger for a: = —1 than a; = 0. As a result, the EOS 
of hypemuclear matter shows similar stiffness for differ- 
ent symmetry energy parameters. Furthermore, the EOS 
from the MDI-Hyp-R interaction is stiffer than that from 
the MDI-Hyp-A interaction, as E particles do not appear 
in the former case as shown in Fig. |6l Our results thus 
show that the hypemuclear matter has a stiffer (softer) 
EOS at lower densities but a softer (stiffer) EOS at mod- 
erate densities with a stiffer (softer) nuclear symmetry 
energy. We note that the results from both x = and 
X ~ —\ are consistent with the constraints on the nu- 
clear equation of state obtained from the analysis of the 
collective flow data in heavy- ion collisions [70j . 



100 




FIG. 9: (Color online) Particle fractions in a hypemuclear 
matter with the presence of hadron-quark phase transition 
from the MDI-Hyp-R interaction with x — ((a) and (c)) 
and X = —1 ((b) and (d)) for the hadron phase and the MIT 
bag model for the quark phase. Results from B^^^ = 180 
MeV ((a) and (b)) and 170 MeV ((c) and (d)) are shown for 
comparison. 

At higher densities in the core of a neutron star, a tran- 
sition from hadron matter to quark matter is expected 
to occur. Here we only consider the hadron-quark phase 
transition for the repulsive EN interaction as it is more 
consistent with the latest empirical information [56l - [60| . 
In Fig. [SI we display the particle fractions of each species 
in the presence of a hadron-quark phase transition, with 
the hadron phase from the MDI-Hyp-R interaction with 
X = and a; = — 1 and the quark phase from the MIT 
bag model with bag constants B^^^ = 180 MeV and 170 
MeV. The particle fractions are weighted by the baryon 
number and corresponding phase fraction. It is seen that 
for B^/* — 180 McV; the phase transition begins at the 
density 0.35 fm~'^ and ends at the density 1.38 fm"-^ for 
a: = 0, and begins at the density 0.26 fm""^ and ends at 
the density 1.38 fm-^ for x = -1, while for B^/^ = 170 
MeV it begins at the density 0.21 fm~^ and ends at the 
density 1.05 fm""^ for a: = 0, and begins at the den- 
sity 0.18 fm^'^ and ends at the density 1.05 fm~'^ for 
a: = — 1. The hadron-quark phase transition thus hap- 
pens at lower baryon number density for a stiffer symme- 
try energy and for a smaller value of B, while the density 
at which the hadron-quark phase transition ends seems to 
depend only on the value of B but not much on the value 
of the symmetry energy parameter x. With a smaller 
value of the bag constant, the hadron-quark phase tran- 
sition both begins and ends earlier. It is also seen that 
d and s quarks occupy a larger fraction than u quarks 
in the mixed phase because of their negative charges, so 
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the lepton fraction decreases while the proton fraction in- 
creases when the hadron-quark phase transition occurs. 
Furthermore, only A hyperons (no other hyperons) ap- 
pear in the mixed phase in our model. The fraction of 
A hyperons is, however, sensitive to the value of the bag 
constant B, and with a smaller value of B its fraction 
becomes smaller. 
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FIG. 10: (Color online) Same as Fig. [9] for the equation of 
state. 

For the EOSs in the presence of the hadron-quark 
phase transition, they are shown in Fig. [TU] with B^l^ = 
180 MeV and 170 MeV. Our resuhs show that the EOSs 
of the mixed phase and the quark phase are both softened 
in comparison with a pure hadron phase. The difference 
between the results from x = and a; = — 1 is not large 
except for the different starting density of the hadron- 
quark phase transition. Since the EOS is more sensitive 
to the value of bag constant than to the value of x when 
the quark degrees of freedom are introduced, the bag 
constant B is thus the main parameter in determining 
the EOS of dense matter as reported in other work [g^] . 
The energy density at the end of the hadron-quark phase 
transition is about 1.2 GeV/fm^ for B^l^ = 180 MeV 
and about 1.6 GeV/fm^ for B^/^ = 170 McV, with the 
former closer to the value of about 1 GeV/fm^ obtained 
from the lattice QCD calculation [Til and extracted from 
heavy- ion collision experiments [72l | . Compared with the 
equation of state constrained by the collective flow data 
in heavy- ion collisions [70| . our results from both values 
of bag constants satisfy the empirical constraint. 



C. Hybrid stars 

In this subsection, we use the MDI-Hyp interaction to 
study the properties of static hybrid stars with spher- 
ically symmetric mass distributions. In particular, we 
calculate the mass-radius relation of a hybrid stars using 



Tolman-Oppenheimer-Volkoff (TOV) equation [TJ, IZ; 



dP 
dr 



r{r — 2mg) 



(30) 



where nig is the gravitational mass inside the radius r of 
the hybrid star given by 



drUg 
dr 



47rr^e(r), 



(31) 



with e{r) being the energy density. Above equations are 
solved by starting from a central energy density e(r = 
0) = Ec and integrating outward until the pressure on 
the surface of the hybrid star defined by r = i? vanishes, 
i.e, P{R) = 0. This gives the radius R, and the total 
gravitational mass of a hybrid star is then given by M = 
mgiR). 

In our calculations, the hybrid star is divided into three 
parts from the center to the surface: the liquid core, 
the inner crust and the outer crust as in our previous 
work [13, In the liquid core it is assumed to be 

the hypernuclear matter or that with the hadron-quark 
phase transition, and the resulting EOSs shown in pre- 
vious subsections are used. In the inner crust, a param- 
eterized EOS of P = a -I- 6e^/'^ is used as in the previous 
treatment [l^, [3]. The outer crust usually consists of 
heavy nuclei and the electron gas, where we use the BPS 
EOS [t^. The transition density pt between the liquid 
core and the inner crust has been consistently determined 
in our previous work [l2l . [l3l | , and for the density which 
distinguishes the edge of inner crust and outer crust, we 
take it to be pout = 2.46 x 10~* fm~'^. The parameters 
a and b can then be determined by the pressures and 
energy densities at pt and Pout- 

In Fig. 1111 we show the mass-radius (M-R) and mass- 
central density (M-pc) relations of hybrid stars from the 
EOS of the hypernuclear matter without the hadron- 
quark phase transition as obtained above based on the 
MDI-Hyp interaction. Results using both attractive and 
repulsive SN interactions are shown, and the results from 
a pure nucleonic approach without hyperons are also 
shown for comparison. The maximum mass obtained 
with hyperons in hybrid stars is seen to be smaller as 
the EOS of the hypernuclear matter is softer than that 
of the nuclear matter. The EOSs at moderate densi- 
ties play an important role in determining the maximum 
mass, as the M-pc relation shows that the neutron star 
mass can not increase even the central density reaches 
a very high value. The mass is larger when a repulsive 
SN interaction is used as a result of the smaller number 
of degrees of freedom and stiffer EOS, as discussed in 
the previous subsection. For smaller values of pc a stiffer 
symmetry energy gives a larger mass as the nuclear mat- 
ter dominates this density range, and the results are thus 
the same as those in our previous work [l^, US]- With 
larger pc, the effect due to the nuclear symmetry energy 
is small, and the mass from a; = is similar to that from 
x = — 1 because of the opposite effects from the symme- 
try pressure and the fraction of hyperons on the equation 
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FIG. 11: (Color online) The hybrid star mass as a function 
of radius (left panel) and central density (right panel) based 
on the MDI-Hyp interaction with a; = and x = —1. Results 
from attractive and repulsive EN interactions are shown for 
comparison. Results from a pure nucleonic approach are also 
displayed. 
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FIG. 12: (Color online) The hybrid star mass as a function 
of radius (left panel) and central density (right panel) in the 
presence of the hadron-quark phase transition. Results from 
the MDI-Hyp-R interaction for the hadron phase with x = 
and X = —1 and the MIT bag model for the quark phase with 
B^^'^ = 180 MeV and 170 MeV are shown for comparison. 
Results from a pure nucleonic approach are also displayed. 



of state. The maximum mass obtained with an attrac- 
tive EN interaction is 1.59Mo for a; = and l.STM© for 
X — —1, where Mq is the solar mass. With a repulsive 
SN interaction the maximum mass increases to 1.65Mq 
for X = Q and 1.67 Mq for x = —1. Our results thus dif- 
fer from those of many previous works [s^, that the 
maximum mass of hybrid stars cannot reach the canoni- 
cal value of IAMq. In our model, the radius of a hybrid 
star with a mass of 1.4AjrQ is 11.2 km for x ~ and 11.9 
km for a; = — 1 for an attractive SN interaction, and 11.9 
km for X ~ and 13.2 km for x ~ ~l for a repulsive SN 
interaction, respectively. 

In Fig. [T^l the M-R and M-pc relations are displayed 
for hybrid stars with the hadron-quark phase transition 
in their liquid core. Again, results from the hadron phase 
with X = and x = —1 and the quark phase with B^^^ = 
180 MeV and B^/^ = 170 MeV are shown, and those from 
a pure nucleonic approach are also shown for comparison. 
The maximum mass for = 180 MeV is I.SOMq for 
x = and 1.46Mq for x = -1, and for B^^-^ = 170 
MeV it is IAGMq for a; = and 1.45M0 for a; = -1, 
respectively. The radius of a standard neutron star with 
a mass of 1.4AfQ is 11.0 km for a- = and 10.8 km for x ~ 
-1 for ^1/4 = 180 MeV, and 10.2 km for a: = and 10.0 
km for X = —1 for S^/'* = 170 MeV. If we further reduce 
the value of B, the hadron-quark phase transition would 
happen at an even lower density, and the maximum mass 
for the hybrid star would correspond to a larger central 
density and a smaller radius. Above results are obtained 
with the Gibbs construction for the hadron-quark phase 



transition. With the Maxwell construction, the radius of 
a maximum mass hybrid star would be larger. Within 
our model, both treatments of the phase transition give, 
however, reasonable masses and radii for hybrid stars. 

Finally, we note that the original MDI interaction does 
not lead to the violation of causality in /3-stable npefi 
matter up to lOpo as shown in Panel (d) of Fig. 2 in 
Ref. [l3|- Since the EOS of neutron star matter is soft- 
ened by the presence of hyperons and quarks, the causal- 
ity condition is still satisfied for the extended MDI inter- 
action used in the present study, and this is confirmed by 
explicit calculations. 

V. SUMMARY 

We have extended the MDI interaction for the nucleon- 
nucleon effective interaction in nuclear medium to include 
the nucleon-hyperon and hyperon-hyperon interactions 
by assuming that they have the same density and mo- 
mentum dependence as that for the nucleon-nucleon in- 
teraction. The parameters were determined by fitting the 
empirical hyperon single-particle potentials in symmetric 
nuclear matter at its saturation density. As an example 
for the application of the extended MDI interaction, we 
investigated the properties of hybrid stars that include 
not only the hyperon degree of freedom but also that of 
quarks by taking into account the phase transition be- 
tween the hadron and quark phases. Our results indicate 
that the extended MDI interaction can give a reasonable 
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description of the hypernuclear matter. We found that 
the EOS of the hypernuclear matter is much softer than 
that of the pure nuclear matter, and that it becomes even 
softer if the hadron-quark phase transition is included. 
The masses and radii of hybrid stars were also studied 
with these EOSs, and they were found to remain rea- 
sonable after including hypcrons and the hadron-quark 
phase transition. We have also studied the effects of at- 
tractive and repulsive EN interactions and different sym- 
metry energies on the hybrid star properties. The results 
show that the appearance of the S hyperon in hybrid 
stars depends sensitively on the sign of the EN interac- 
tion with a repulsive EN interaction giving a higher criti- 
cal density for the appearance of E hyperons. In addition, 
a stiffer symmetry energy usually leads to a larger frac- 
tion of hyperons in the hypernuclear matter. We further 
found that both the low-density boundary of the hadron- 
quark phase transition and the EOS at high densities in 
hybrid stars are more sensitive to the bag constant than 
to the stiffness of the nuclear symmetry energy at high 
densities. This extended MDI interaction, which gives 
an isospin- and momentum-dependent in-medium effec- 
tive interactions for the baryon octet, will also be very 



useful in transport models that simulate heavy-ion reac- 
tions in future radioactive beam facilities, particularly at 
the FAIR/GSI energies. 
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